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Iron abundance in H ii regions 
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Abstract. Optical CCD spectra are used to determine the Fe abundances at several positions inside seven 
bright Galactic Hll regions. The observed [Felll] line ratios are compared with the predictions of different sets 
of collision strengths and transition probabilities for this ion to select the atomic data providing the best fit to 
the observations. The values found for the Fe^""" and Fe""" abundances, along with ionization correction factors 
for the contribution of Fe^+, obtained from available grids of photoionized models, imply that the Fe/O ratio in 
the ionized gas is between 2% and 30% of solar. The Fe abundances derived for each area are correlated both 
with the degree of ionization and the colour excess. A possible explanation is suggested, namely the presence of 
a population of small grains, probably originating from the fragmentation of larger grains. These small grains 
would release Fe atoms into the gas after the absorption of energetic photons; the small grains surviving this 
destruction process would be swept out of the ionized region by the action of radiation pressure or stellar winds. 
An indication of a further and more efficient destruction agent is given by the high Fe abundance derived for a 
position sampling the optical jet H 399 in M20, where dust destruction due to shock waves has presumably taken 
place. 
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1. Introduction 

The first reliable measurement of the Fe abundance in an 
Hii re gion was carried out for M42 by Osterbrock et al. 
( |1992D , who found Fe/H to be underabundant by a factor 
] of 12 with respect to the solar value, confirming the previ- 
■ ous estimate by Olthof & Pottasch (1975). Similar results 
' were obtained for M17 and M8, for which Pcimbert et 
al. ( 1993| ) found underabundance factors of 4 and 20, re- 
spectively. These values are consistent with the expected 
depletion of Fe on to the dust grains known to coexist 
with the ionized gas (Miinch & Persson 1971; Osterbrock 



1989), but the derived Fe abundances are significantly 



above the typical values for the undisturbed dense inter- 
stellar medium, ~ 100 below solar (Jenkins 1987] ). This 
difference in the Fe depletion factors in H ii regions and 
the interstellar medium is not readily explained, since re- 
fractory dust grains are expected to survive inside Hii 
regions (Osterbrock p^ ). 

A possible explanation was offered in Rodriguez 
( |1996| ), where the gaseous Fe abundances at different po- 
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sitions within seven Galactic Hii regions were shown to 
be correlated with the degree of ionization, and this cor- 
relation was interpreted by relating the energetic photons 
with the release of Fe atoms from dust grains. This paper 
presents a more detailed and complete analysis of the data 
used in Rodriguez ( 1996 ) to determine the Fe abundances 
in M42, M43, M8, M16, M17, M20 and NGC 7635. The 
present analysis relies on recent atomic data for the de- 
termination of physical conditions and ionic abundances, 
and uses a somewhat different approach to derive the Fe 
abundance, leading to abundance values lower by about 
40% than those presented in Rodriguez (1996). The corre- 
lation of the Fe abundance with the degree of ionization is 
confirmed to hold, as well as the further correlation with 
the colour excess, also previously presented in Rodriguez 
( 19961 ). The implications that both results have for dust 
in H II regions are discussed. 



2. The data 

The results in this paper are based on the line in- 
tensities measured by Rodriguez (|199|, [l999b| ) in the 
AA4200-8800 spectra of several positions within seven 
Galactic Hii regions: M42, M43, M8, M16, M17, M20 
and NGC 7635. This work centres on the Fe lines, 



but the values derived by Rodriguez ( |l999bD for some 
physical conditions and ionic abundances are also used: 
iVe[Sii], the densities implied by the intensity ratio of 
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[Sii] lines /(A6716)//(A6731); To[Nii], the temperatures 
derived from the intensity ratio of [N ii] hnes /(A6548 + 
A6583)/J(A5754); and the ionic abundance ratios N+/H+, 
0+/H+, 0++/H+ and 0/H = 0+/H+ + 0++/H+. 

The ionization potentials for Fe°, Fe+, Fe++ and Fe'^+ 
are 7.9, 16.2, 30.6 and 54.8 eV, respectively. Considering 
the excitation characteristics of the objects in the sample 
and the absence of He ii A4686 emission from all the neb- 
ular spectra (the ionization potential for He+ is 54.4 eV) 
Fe^+ is not expected to be further ionized, and only Fe+, 
Fe"*"^ and Fe^+ should have appreciable concentrations 
within the ionized gas. Several [Feiii] and [Fell] lines can 
be measured in the optical spectra of Hii regions, but 
no optical [Feiv] lines have been detected in Hii regions 
and only one line has been measured in the UV spe ctra o f 
an Hii region: [Feiv] A2837 in M42 (Rubin et al. |l997|) . 
Therefore, the values derived here for the Fe abundance 
are based on the Fe+ and Fe++ abundances and ionization 
correction factors for the contribution of Fe'^+ to the total 
abundance derived from available photoionization models. 

3. [Feiii] lines 

Up to twelve [Feiii] lines, previously listed and measured 
in M42 by Osterbrock et al. ( |1992D and belonging to (IF), 
(2F) and (3F) multiplets, were detected in each object. At 
the achieved spectral resolution of 4 A, [Fe ill] A4986 and 
[Felli] A4987 are blended; [Feiii] A4607 is co ntami nated 
by On A4609 and Nii A4607 (Esteban et al. |1998D ; the 
intensity measured for [Feiii] A4778 might have a small 
contribution from Nil A4780 (Esteban et al. 1999a), and 
[Feiii] A4658 and [Feiii] A5270 are appreciably blended 
in some objects with On A4661 and [Feii] A5273, respec- 
tively. On the basis of these considerations, [Feiii] A4607 
has been excluded from the abundance determination, 
[Feiii] A4778 has been considered only in those positions 
where the abundance implied by this feature is consistent 
with the mean abundance implied by the other lines, and 
[Feiii] A4658 and [Feiii] A5270 have been used only for 
those objects where the contamination of these lines is 
negligible. 

The reddening-corrected intensity ratios of the [Fe ill] 
lines relative to H/3 are listed in Ta ble 1 for all the posi- 
tions studied (see Rodriguez 1999b for the values of the 
parameters used in the extinction correction) . The extinc- 
tion corrections applied to the relative [Feiil] line inten- 
sities are very small, a few per cent in most cases, with a 
maximum of 30% for [Felll] A5270 in M17. The accuracy 
of the line intensities in Table 1 may vary from 10-20% 
for the positions in M42 to more than 50% for some cases, 
in particular for those positions where only one [Fe ill] line 
could be measured. 

4. Atomic data for [Fein] 

Keenan et al. (|1992|) derived the relative populations of 



transition probabilities of Garstang (1957). Further calcu- 
lations of the collision strengths for this ion are reported 
by Zhang (1996), whereas new values for the transition 
probabilities for the forbidden transitions have been de- 
rived by Nahar & Pradhan ( |1996| ) and Quinet ( |1996D . The 
new effective collision strengths, when compared with the 
previous values of Berrington et al. (1991), are enhanced 
by a factor of up to 2, whereas the three available sets of 
transition probabilities differ from each other by ^ 30% 
for most transitions, with differences up to an order of 
magnitude or more for some transitions. 

In order to compare the results implied by the most re- 
cent atomic data with the previous calculations of Keenan 
et al. ( 1992 ), the level populations for Fe++ have been 
calcul ated with the collision strengths reported by Zhang 
( |1996|) and the two recent sets of tran sition probabilities 
(Nahar & Pradhan 1996 ; Quinet 1996). All coUisional and 
downward radiative transitions for the 34 levels of the 3dP 
configuration have been considered; test calculations re- 
stri cted t o the lowest 17 levels considered by Keenan et 
al. (|1992D are found to agree to within 10% with the 34- 
level calculations. 

In Fig. H some observed and predicted [Fe ill] line ra- 
tios are presented as a function of A'c [S ii] (Rodriguez 
1999b). The values found for To [Nil] cover the range 



rc[Nii] = 7800-10 800 K (Rodriguez |l999b|), but the pre- 
dicted line ratios for any temperature in this range are sim- 
ilar, and only the calculations for Tc = 10 000 K are shown 
in Fig, p]. T he [Feiii] line ratios measured by Osterbrock 
et al. ( |1992[ ) and Esteban et al. ( |1998| , |l999aD in M42 and 
M8 are also presented in this figure; note that the mea- 
surements of Esteban et al. ( [1998 , 1999a), performed at 



the lowest 17 levels of the 3dP configuration in Fe "'"^ using 
the collision strengths of Berrington et al. (1991) and the 



high-resolution, can be considered the most reliable. 

t igures |l|a, |l|b, |l|d and 0e show some of the line ra- 
tios for which the different atomic parameters lead to the 
greatest differences. In general, all the atomic data fit the 
measured line ratios reasonably well, but the results in 
Figs. ^ and |^ suggest that t he cal culations based on 
the collision strengths o f Zha ng ( 1996 ) and the transition 
probabilities of Quinet ( 1996 ) lead to a better fit. 

The best way to judge the accuracy of the atomic data 
is to consider the dispersion in the abundances implied by 
the different lines. The most accurate [Fe ill] intensities are 
probably those measured by Esteban et al. ( 1998| ) in M42. 
The abundances and dispersions implied by the intensities 
measured in their position 1 for the ten lines in Table 1 and 
the three combinations of atomic data considered here are 
the fo llowing: Fe++/H+ = (3.5± 0.8) x 10"'^ (Keenan et al. 
19921 ), (2.5 ± 0.3) X 10"'^ (Zhang p^; Na har & Pradhan 
m§ and (2.5 ± 0.2) x 10"^ (Zhang |1996| ; Quinet |l99|). 
When the [Fe III] intensities measured in other regions are 
considered, similar results are obtained: the oldest atomic 
data lead in most cases to higher dispersions and 10 to 50% 
higher mean abundance values. The most recent collision 
strengths lead to essentially equal mean abundance values 
independently of the transition probabilities used, and the 
dispersions are similar for the two recent sets of transi- 
tion probabilities, although the transition probabilities of 



M. Rodriguez: Iron abundance in Hll regions 
Table 1. [Feiii] line intensities and Fc++ abmidance 
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Object /(A)/7(H/3) x 100=^ Fe++/H+ ± a 







A4658 


A4702 


A4734 


A4755 


A4769 


A4778 


A4881 




A5270 


A5412 








M42 A-1 




0.77 


0.255 


0.096 


0.142 


0.086 


0.081: 


0.32 


0.056 


blend 


0.029 


(3.3±0.4] 


xlO" 


7 


A yr /I o A n 
M4z A—Z 




0.69 


0.208 


0.082 


0.132 


0.071 


0.061 


0.263 




blend 


0.026 


(2.2±0.4 


X 10~ 


7 


M42 A-3 




0.64 


0.189 


0.065 


0.113 


0.065 


0.029 


0.255 




blend 


0.015 


(2.5±0.4 


X 10" 


7 


M42 A-4 




0.72 


0.23 


0.090 


0.153 


0.074 


0.064: 


0.289 




blend 


0.021 


(2.7±0.4 


X 10~ 


7 


ivr42 A-fi 




0.72 


0.20 


0.074 


0.132 


0.047 


0.042 


0.25 




blend 


0.015: 


(3.6±0.6 


X 10~ 


7 


M42 A-fi 




0.685 


0.21 


0.06 


0.14 






0.249 




blend 




(2.5±0.2^ 


X 10~ 


7 


M4z _D-i 




0.87 


0.251 


0.088 


0.179 


0.081 


0.052 


0.352 


0.065 


blend 


0.031 


(4.2±0.7 


X 10~ 


7 


M42 R-2 




0.905 


0.248 


0.110 


0.186 


0.094 


0.047 


0.364 


0.05: 


blend 


0.023 


(4.7±0.8" 


X 10~ 


7 


M42 R-'H 




0.85 


0.219 


0.067 


0.18 


0.074 


0.052 


0.309 


0.053: 


blend 


0.018 


(4.2±1.0 


X 10" 


7 


M42 B-4 




1.34 


0.384 


0.149 


0.262 


0.128 


0.083 


0.518 


0.118 


blend 


0.053 


(6.7±1.0 


X 10" 


7 


M42 B-5 




0.894 


0.25 


0.07 


0.14 


0.08 


0.05 


0.287 


0.11 


blend 




(4.8±1.1 


X 10" 


7 


M42 B-6 




0.90 


0.26 


0.08 


0.146 


0.10 


0.06 


0.33 


0.08 


blend 


0.038 


(4.9±0.7 


X 10" 


7 


M43-1 




0.8 


0.16 




0.19 






0.24 


0.178: 


blend 




(8.6±1.9 


X 10" 


7 


M43-2 




0.81 


0.23 




0.33; 






0.249 


0.18: 


blend 




(8.6±0.9 


X 10" 


7 


M43-3 




0.6 






0.12 








0.21 


blend 




(5.5±0.9) 


xlO" 


7 


M43-4 




0.83 


0.21 




0.19 






0.209 


0.164: 


blend 




(9.4±1.2] 


xlO" 


7 


M43-5 




0.94 


0.18 




0.19 






0.32 


0.17: 


blend 




(1.0±0.2] 


XlO" 


6 


M8-1 




0.57 


0.18 


0.07 


0.11 


0.06 




0.25 


0.07 


blend 




(3.2±0.4] 


XlO" 


7 


M8-2 




0.57 


0.18 


0.05 


0.065 






0.21 


0.06 


blend 




(3.1±0.9) 


XlO" 


7 


M8-3 




0.64 


0.19 




0.13 


0.04 




0.23 




blend 




(3.7±0.7] 


XlO" 


7 


M8-4 




0.55 


0.19 


0.05 


0.10 


0.037 


0.02 


0.178 


0.052 


blend 




(3.2±0.8) 


XlO" 


7 


iViO rJ 






n 9'^- 

u.zo. 




nil 






n 1 ^ 


n no 


blend 




1 r> c 1 1 I ^ 


V 1 n~ 


7 


M8-6 




0.57 


0.11 










0.11 


0.10 


blend 




(2.3±0.5) 


XlO" 


7 


M16-1 




0.17 






0.06 






0.044 




blend 




(1.5±0.8) 


XlO" 


7 


Ml 6-2 


















0.31 






2.0xl0~" 






M17-1 




blend 












0.039 




0.125 




(1.2±0.5) 


XlO" 


7 


Ml 7-2 




blend 






0.08 






0.07 




0.12 




(1.5±0.5] 


XlO" 


7 


Ml 7-3 




blend 












0.07 




0.095 




(1.3±0.3; 


xlO" 


7 


M20-1 




1.0 


0.18 




0.19 






0.30 


0.29: 


blend 




(7.8±1.9] 


XlO" 


7 


M20-2 




0.3 












0.04 


0.23 


blend 




(2.2±0.6) 


XlO" 


7 


M20-3 




0.4 




















2.9x10-" 






NGC 7635- 


-1 


0.14 


0.06 










0.122 


0.04 


blend 




(1.7±0.6] 


XlO" 


7 


NGC 7635- 


-2 


0.39 


0.11 




0.08 






0.18 


0.046 


blend 




(2.5±0.2] 


XlO" 


7 


NGC 7635- 


-3 


0.38 






0.07 






0.11 


0.04 


blend 




(2.1±0.3) 


XlO" 


7 



[Feiii] A4658 and [Feiii] A5271 are blended with On A4661 and [Feii] A5273, respectively. The line intensities marked 
with colons imply abundance values that deviate by more than 50% from the mean value derived from the other lines 
and have been excluded from the abundance determination 



Quinet (1996) lead in general to somewhat lower disper- 
sions. Therefore, the collision strengths of Zhan g ( 1996 ) 
and the transition probabilities of Quinet (1996) will be 
used hereafter to derive the Fe^"'' abundances. 



5. Fe++/H+ 

The last column in Table 1 shows the values of Fe++/H+ 
derived at each position from the measured [Fe ill] i nten- 
sities by using the collision strengt hs of Zhang ( 1996| ), the 
transition probabilities of Quinet ( |l996|) , the Hi emissiv- 
ities of Hummer & Storey ( 1987| ) and th e physical condi- 
tions Te[Nii], iVe[Sii] (Rodriguez |l999bD . The dispersions 



in the abundances implied by the different lines considered 
are also indicated in Table 1; it can be seen that most of 
them are quite low, in the range 10-20%. 

The values derived forFe^+ZH"*" are lower than those 
presented in Rodriguez ( |1996[ ) by 10-70%. Up to 20% in 
these differences is due to the new values of the physi- 
cal conditions implied by updated atomic data (Rodriguez 



1999b); the remaining differences arise from the use of the 



most recent atomic data for Fe"*""*". 
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Fig. 1. Several observed [Feiii] line ratios are presented as a function of density, A'o[Sii], and compared with the 
collisional excitation predictions for To = 10 000 K. The atomic data use d in t he calculations are id entifie d in panel e, 
where Z96, Q96, NP96, BZ9 1 and G57 stand for Zhang (|l996|) , Quinet (|l996| ), Nahar & Pradhan (|l99^) , Berrington 
et al. ( 1991 ) and Garstang ( 1957 ), respectively. The symbols representing the different Hii regio ns are identified in 
panel f, where OTV 92, EP98 and EP99 identify the measurements performed by Osterbrock et al. (1992) and Esteban 
et al. (|l998i |l999a|) 



6. [Fell] emission and Fe+ abundance 

Many lines of [Fell] have been identified in the optical 
spectra of Hii regions (see, for example, Osterbrock et 
al. |1992| ; Rodriguez |1996| ; Esteban et al. |1998| , |l999aD . 
Most of these lines are s everely affected by fluorescence 
effects (Rodriguez 1999a ), but an estimate of the con- 
tribution of Fe"*" to the total Fe abundance can be ob- 
tained from the intensity measured for [Fell] A8617, which 
is alm ost insensitiv e to th e effects of UV pumping (Lucy 



1995| ; Baldwin et al. |1996[ ) . The reddening-corrected inten- 



sities of [Fell] A8617 relative to Pal2 or Pal3 are listed in 
Table 2; the applied extinction corrections are below 10%, 
and the uncertainties in the relative intensities may vary 
widely from ~15% in M42 to more than 50% in M20-2, 
M16 and M17. The values of Fe"'"/H+ have been derived 
from the intensities in Table 2 by us ing the emissivities 
for [Fell] (Baut ista & Pradhan 1996D and Hi (Hummer 
& Storey 1987) interpolated to the physical conditions 
To [Nil] and iVc[Sii]. 

In general, the values of Fc+/H+ are much lower than 
those found for Fe+^/H+j as can be expected from the 
low ionization potential of Fe"*" (16.2 eV): most of the 
positions show Fe"'"/Fe+"'" ^ 0.3 (see Table 2). One ex- 
ception is M43-3, where the Fe++/H+ abundance ratio is 



below those derived for the other positions in M43, and 
the value of Fc++/H+ + Fc+/H+ turns out to be very sim- 
ilar. The most notable exceptions are M20-1 and M20-2, 
where Fe+/Fe++ ~ 1. The spectra obtained for M20 are 
very noisy in the wavelength range around [Fell] A8617 
and the intensity measured for this line is uncertain by a 
factor of 2 in M20--2, but in position M2G-1 aU the [Fell] 
and [Fe ill] lines are relatively strong and the intensity of 
[Fe ii] A8617 has a much lower uncertainty. Position M20-1 
samples the optical jet HH 399, where the conditions could 
differ from those found in other regions (see the discussion 
in Sect. 7.1). 



It has been suggested (Bautista & Pradhan 1998, and 
references therein) that most of the [Fe ii] emission in M42 
and other Hll regions originates in partially ionized lay- 
ers with densities N^, ~ 10^ cm~'^, much higher than those 
measured in the fully ionized zones, where iVe < 10^ cm~^, 
but several arguments have been presented against this 



hypothesis (Baldwin et al. 199(: ; Esteban et al. 1999b 



Rodriguez 1999a; Verner et al 



2000). In any case, the 



[Fe li] lines measured in the near- infrared spectrum of M42 
have been foun d to arise in regions of m oderate density 
(Marconi et al. 1998 ; Luhman et al. 199^ ), similar to the 
densities measured from the usual diagnostics based on 
[Sii], [On] or [CI III] lines, and the [Fell] A8617 intensity 
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Table 2. [Fell] A8617 intensities and Fe+ abundance 



Object 




J//(Hi) 


Hi 


Fe+/H+ 


Fe+/Fe++ 


M42 A-1 




0.094 


Pal2 


7.8x 10~* 


0.24 


ivr49 A -9 




0.09 


Pal2 


5.9x 10~* 


0.27 


M42 A-3 




0.030 


Pal2 


2.4x10"* 


0.10 


M42 A-4 




0.061 


Pal2 


4.5x10"* 


0.17 


M42 A-5 




0.048 


Pal2 


4.5x10"* 


0.13 


M42 A-6 




0.038 


Pal2 


2.9x10"* 


0.12 


M49 R-1 




0.088 


Pal2 


7 gx 10"* 


0.19 


M42 B-2 




0.044 


Pal2 


4.1x10"* 


0.09 


M42 B-3 




0.055 


Pal2 


5.2x10"* 


0.12 


M42 B-4 




0.159 


Pal2 


1.4x10"^ 


0.21 


M42 B-5 




0.08 


Pal2 


8.5x10"* 


0.18 


M42 B-6 




0.101 


Pal2 


9.8x10"* 


0.20 


M43-1 




0.095 


Pal2 


1.7x 10"^ 


0.20 


M43-2 




0.157 


Pal2 


2.7x10"^^ 


0.31 


M43-3 




0.21 


Pal2 


3.6x10"^ 


0.65 


M43-4 




0.085 


Pal2 


1.6x10"^ 


0.17 


M43-5 






Pal2 


1.4x 10"^ 


0.14 


M8-1 




0.021 


Pal2 


2.0x10"* 


0.06 


M8-2 




0.027 


Pal2 


3.0x10"* 


0.10 


iVio o 






Pnl 9 

jra_L^ 


9 c: V 1 


U.U ( 


M8-4 




0.013 


Pal2 


1.5x10"* 


0.05 


M8-5 




0.019 


Pal2 


2.6x10"* 


0.07 


M8-6 




0.02 


Pal2 


2.2x10"* 


0.10 


M16-1 




0.015 


Pal2 


1.9x10"* 


0.13 


M16-2 












M17-1 




0.009 


Pal2 


l.lx 10"* 


0.09 


Ml 7-2 












Ml 7-3 




0.016 


Pal3 


1.5x10"* 


0.12 


M20-1 




0.5 


Pal2 


7.7x10"'^ 


0.99 


M20-2 




0.1 


Pal2 


1.9x10"^ 


0.86 


M20-3 












NGC 7635- 


-1 


0.051 


Pal2 


5.9x10"* 


0.35 


NGC 7635- 


-2 


0.066 


Pal2 


6.7x10"* 


0.27 


NGC 7635- 


-3 


0.06 


Pal3 


5.1x10"* 


0.24 



measured in M42, when compared with the intensities of 
these infrared hues, impHes a common origin at moderate 
densi ties (see, for example, Fig. 5d in Bautista & Pradhan 
Hence, the calculation of the Fe^/H"*" abundance 



1998 



ratio from [Fell] A8617 for the physical conditions To [Nil] 
and Ne [S ii] seems quite reasonable. 

Another source of uncertainty in the values of Fc+/H+ 
in Table 2 arises from the atomic data used for Fe+ - par- 
ticularly the collision strengths. The Fe+ abundances have 
been calculated with the emissivities derived by Bautista 
& Pradhan ( |1996[ ) using their own collision strengths, al- 
though Bautista & Pradhan ( 199^ ) consider more accu- 
rate the collision stren gths o f Pradhan & Zhang ( 1993 ) 
and Zhang & Pradhan ( 1995 ). However, this other set of 
collision strengths would lead to very similar values (~30% 
lower) for the Fe+ abundance derived from [Fell] A8617. 
The collision strengths presented by Pradhan & Zhang 



(1993) and Zhang & Pradhan (1995) can also be affected 
by uncertainties (Rodriguez 1999a; Oliva et al. 1999), but 
the contribution of Fc+ to the total abundance is low and 
these uncertainties will not affect significantly the derived 
Fe abundances. The Fc abundance could be alternatively 
derived from the values of Fe++/H+ and ionization cor- 
rection factors. This is the approach followed in Rodriguez 
( 1996D , and leads to very similar results. 



7. Fe abundance 

The values of log (Fc++/H+ + Fc+/H+) are presented as 
functions of the degree of ionization, log (0+/0++) and 
log (N+/N), in Fig. |, where log (N/H) = -4.35 has been 



assumed to hold for all the nebulae (Rodriguez 1999b). 
Figure || shows that Fe+"'"/H+ -I- Fe"'"/H+ decreases with 
increasing degree of ionization for the positions in M42 
and M17, reflecting the expected increment in the con- 
tribution of Fe^+ to the total abundance. However, the 
regions with lower degree of ionization - excluding M43 
- show unexpectedly low values of Fc++/H+ + Fc+/H+, 
that cannot be due to an error in the derived Fe+/H+ 
abundance ratio, as argued above and in Sect. 7.1. 

Figures ||a and ||b also show an indication of the ex- 
pected - i.e. model predicted - behaviour of the ion- 
ization fractions x(Fe+"'") -I- a::(Fe"'") as a function of 
a;(0+)/x(0++) and a;(N+), scaled to a total abundance 
log(Fe/H) = —5.8. Two grids of ionization models have 
been used: those from Stasihska (1990) for solar metallic- 
ity Zq and for Zq/2, and those from Gruenwald & Viegas 
( |1992| ) for Zq and Zq/3 (see Rodriguez p^99b| for further 
details on the selected models). The vertical dotted lines 
in Figs. I^a and show the range covered by the individ- 
ual photoionization models; a further dotted line shows 
the assumed ionization correction factors, also scaled in 
these figures to log(Fe/H) = —5.8, and given by: 



Fc/H = 1.1 X (0/0+) X (FC++/H+ + Fc+/H+) , (1) 
Fe/H = 0.94 x (N/N+) x (Fc++/H+ + Fc+/H+) . (2) 

These relations provide the correction for the contri- 
bution of Fe'^"'' to the Fe abundance in terms of the ionic 
and total abundances of O and N. According to the mod- 
els, Eqs. (|^) and have accuracies better than ±0.2 and 
±0.1 dex, respectively. Fe^+, O^"*" and N++ are formed at 
30.6, 35.1 and 29.6 eV, and none of these ions is expected 
to be further ionized. Hence, the ionization correction fac- 
tor based on N should be more reliable, but whereas the 
total O abundance can be derived from the 0+ and 0++ 
abundances using optical observations, the total N abun- 
dance has to be estimated from N+. Therefore, the results 
for both Fe/N and Fe/0 are presented here, and Fe/H is 
obtained from Fe/0 x 0/H. 

The real accuracy of Eqs.(0) and (H) is difficult to es- 
timate. Even if the ionization and recombination cross- 
sections used by the models are reliable, the predicted 
ionic concentrations may show systematic departures from 
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Fig. 2. The mean values of Fe+/H+ + Fe++/H+ are pre- 
sented as functions of the degree of ionization given by 
0+/0++ (a) and N+/N (b), with log(N/H) = -4.35. 
The relations given by Eqs.(|l|) and @ are also shown, 
scaled to log(Fe/H) — —5.8; The vertical dotted lines 
show the range covered by the individual models (refer- 
ences in the text). It has been assumed that for those 
positions where [Fell] A8617 could not be measured (see 
Table 2), Fe+/II+ has a negligible contribution 



those prevailing in the observed objects due to the neces- 
sarily unrealistic radiation fields and nebular structures 
assumed in the modeling. The dependence of the results 
on the radiation field and the nebular structure can be 
minimized by comparing, for observations and models, 
ions whose ionization potentials are similar, as the ions 
of O and N are compared with the Fc ions in Eqs.(|l]) and 
(^. The dispersion shown by the models around these 
relations in Eqs.(^ and (^, may then reflect real second- 
order variations in the spectral distribution of the radi- 
ation field, the gas density or its metallicity. The con- 
centrations of those ions with low ionization potentials 
(such as 0°, Fe"*" or even 01+ and He+) or relatively high 
ones (such as Ar^+) are, however, severely affected by the 
unrealistic radiation fields or nebular structures assumed 
by models, but the relative ionization fractions of those 
ions whose concentrations depend on the flux density of 
photons with energies of 30-35 eV - as in Eqs.(|l|) and 



(||) ~- seem quite reliable (Rodriguez 1999b). Furthermore, 
the Fe abundances will be more reliable when obtained 
from Fc+/H+ + Fc++/H+ than from Fe++/H+. The re- 
lations implied by Eqs.(|l]) and are also similar to the 
results obtained from individual ionization models. The 



photoionization models presented by Baldwin et al. ( 1991 ) 
and Rubin et al. ( |199l4 |1991bD for two regions in M42 
agree with Eqs.(j^) and (^ to within 0.25 dex. Finally, the 
model derived by Bautista & Pradhan (1998) also for M42, 
which includes the most recent values for the photoioniza- 
tion cross-sections and recombination-rate coefficients for 
the Fe ions, agrees more closely with Eq.(0). In any case, 
to confirm the trends in Eqs.(0) and (§), further calcu- 
lations involving these new ionization and recombination 
cross-sections for different model assumptions would be of 
interest. 

The values derived for Fe/N, Fe/0 and Fe/H are 
listed in Table 3 and shown as a function of degree of ion- 
ization in Fig. 1^. The results are compared with the solar 
abundances in the right hand axes of Fig. |^, where the so- 
lar values are log(Fe/N)o = -0.47, log(F e/0)Q = -1.37 
and log(Fe/H)Q = -4.50 (Grevesse et al. |1996| ). The de- 
pletion factors implied by the Fe/N and Fe/0 abundance 
ratios are broadly consistent, especially if it is taken into 
account that whereas no significant amount of N is ex- 
pected to be present in dust (Sofia et al. 1994), about 
25% of the O abundance (^ 0.1 dex) may be depleted 



in refractory dust grains (Cardelli et al. 1996). However, 
the depletion factors implied by Fe/H are much lower, 
probably reflecting the difficulties related to the choice of 
the correct reference abundances. Several recent studies 
of the abundances of different elements in the interstellar 
medium and in B stars suggest that the abundances of 
O, C, N and Kr rela tive t o H are 2/3 th e sola r val- 



ues (e g. Sn ow & Witt |1996| ; Cardelli & Meyer |1997t Sofia 
et al. 1997), and similar underabundances for O, N and 



other elements are found in Rodriguez (1999b) for all the 
objects studied here. If this result holds also for the Fe/H 
abundance ratio, the most reliable depletion factors must 
be those implied by the Fe/N or Fe/0 abundance ratios. 

The Fe abundances presented here are in general about 
0.2 dex lower than the values derived by other authors (e.g. 
Osterbrock et al. |1992| for M 42; Es teban et al. |1998| , |l999a 
for M42 and M8; Rodriguez |1996D , reflecting the different 
atomic data and the slightly different ionization correc- 
tion factors used. As an example, if the line intensities 



measured by Osterbrock et al. (1992) in M42 are anal- 
ysed with the same atomic data used here and the same 
procedure followed, their derived log(Fe/H) = —5.57 (for 
0+/0++ ~ 1) would change to log(Fe/H) = -5.83 if 
their ionization correction factor, Fe^^/Fe++ = 2, is used, 
and to log(Fe/H) = —5.92 if Eq.(P is used instead. 

The values of the Fe abundance listed in Table 3 are 
highly sensitive to the [Felll] atomic data (in particular 
the collision strengths) a nd to the assumed ionization cor- 
rection factors. Zhang ( 1996 ) estimates uncertainties of 
10-20% for the collision strengths used here to derive the 
Fe^"*" abundances so that, even if the uncertainties are 
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Fig. 3. The values of Fe/N (a), Fe/0 (b) and Fe/H (c), 
as derived from Eq.([|), Eq.(^ and the relation Fe/H = 
Fe/0 X 0/H. The values of Fe/H implied by the Fe/N 
abundance ratio can be obtained from: log(Fe/H) = 



log(Fe/N) - 4.35 (Rodriguez |l999bD 



somewhat higher, the contribution of Fe'^'*' to the total 
abundance seems to be the main source of uncertainty. 

New atomic data for the ionization and recombina- 
tion rates of the Fe ions, or for the rates of their charge- 
exchange reactions, could change the predicted variation 
of x(Fe+) + a;(Fe++) with a;(0+)/a:{0++) or a;(N+), but 
the effects can only be significant for the higher ioniza- 
tion regions where the contribution of Fe'^"'" to the total 
abundance dominates. For those regions, namely the po- 
sitions in M42 A and M17, it might not be surprising to 



find Fc abundances similar to those in M42 B, so that the 
trend in Fig. |^ would level off for 0+/0"'""'" < 1. However, 
the trend followed by the Fe abundances in the regions 
of lower ionization is established by the values derived 
for Fe++/H+ -I- Fe''"/H+ (see Fig. |^) and seems quite ro- 
bust. As discussed above, the real values of the depletion 
factors implied by the derived abundances, are somewhat 
uncertain, but nevertheless they imply that most of the 
Fe atoms remain in the dust. The different degrees of 
depletion factors found for the different objects seem to 
be well established, and the dependence of the results on 
the ionization degree offers an explanation, since it relates 
those photons with energies <; 35 eV to the release of Fe 
atoms from dust|^. The implications of this result will be 
discussed in Sect. 8. 

7.1. Comments on M43 and M20-1 

The Fe abundance ratios derived for M43 and M20-1 show 
significant deviations from the relationship with the de- 
gree of ionization followed by the other objects. M43 is 
the only region in the sample which is excited by a B- 
type star (NU Ori), its spectrum is greatly affected by 
dust-scattered stellar light, and is located several arc min- 
utes to the north of M42, belonging to the same molecular 
complex, although it is an independent H ii region (Goudis 



1982). The Fc/H abundance ratio in M43 seems reliable 
since it is given by the calculated values for Fe+/H+ 
and Fe'^'^/H'*' and does not depend on ionization correc- 
tion factors. A contribution of fluorescence effects, which 



are very important in M43 (Rodriguez 1999a), to the 
[Fell] A8617 intensity cannot be an explanation, since 
when the Fe abundances are derived just from [Fc ill] emis- 
sion, which is not expected to suffer from fluorescence ef- 
fects (Lucy 1995; Bautista & Pradhan 1998D , the same de- 
viation is found for M43 (Rodriguez 1996 ). The relatively 
high Fe abundance in this object is not readily explained; 
maybe it is just an indication of the dependence of the Fe 
depletion on other factors apart from the radiation field: 
age, density, previous history, etc. 

The deviation of M20-1 from the relationships in Fig. |^ 
could be precisely reflecting the action of another factor. 
Position M20-1 is located at the tip of HH 399, as men- 
tioned above, a jet extending into the ionized gas from a 
molecular column in M20 (Cernicharo et al. 199S; Hester 
et al. 1999). The jet is fully ionized and its spectrum is very 
similar to those measured for M20-2 and M20-3, yielding 
similar values for the physical conditions, chemical abun- 
dances and ionization fractions for all elements excluding 



^ The degree of ionization is proportional to th e number of 
energetic photons per electron (Osterbrock 1989), and since 



the distributions of gas and dust within a nebula are expected 
to follow a roughly similar pattern (see Miinch and Persson 



1971 for M42), a higher gas density is likely to be associated 
with a higher dust concentration. Thus, the correlation can be 
interpreted as a relationship between the number of energetic 
photons available per dust grain and the amount of Fe released 
from dust. 
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Table 3. Fc abundances 



Object 




Fc/N 




Fc/O 




Fe/H=' 




M42 A-1 




2.4x10" 


2 


3.7x10" 


3 


1.0x10" 


6 


lVf4-9 A -9 




1.7x 10" 


2 


3.5x 10" 


3 


8. Ox 10" 


7 


M42 A-3 




4.6x10" 


2 


6.3x10" 


3 


1.8x10" 


6 


M42 A-4 




2.7x10" 


2 


4.5x10" 


3 


1.1x10" 


6 


M42 A-5 




3.8x10" 


2 


5.5x10" 


3 


1.5x10" 


6 


M42 A-6 




4.0x10" 


2 


9.0x10" 


3 


1.8x10" 


6 


M42 B-1 




2.2x 10" 


-2 


3.2x10" 


3 


9.4x10" 


7 


M42 B-2 




2.8x10" 


2 


3.7x10" 


3 


1.0x10" 


6 


M42 B-3 




2.3x10" 


2 


3.1x10" 


3 


9.2x10" 


7 


M42 B-4 




2.7x10" 


2 


4.1x10" 


3 


1.2x10" 


6 


M42 B-5 




2.7x10" 


2 


4.1x10" 


3 


1.0x10" 


6 


M42 B-6 




2.4x10" 


2 


3.2x10" 


3 


1.1x10" 


6 


M43-1 




2.3x 10" 


2 


4 5x 10" 


3 


1.2x 10" 


6 


M43-2 




2.6x10" 


2 


4.6x10" 


3 


1.3x10" 


6 


M43-3 




2.1x10" 


2 


2.6x10" 


3 


1.1x10" 


6 


M43-4 




2.3x10" 


2 


4.2x10" 


3 


1.2x10" 


6 


M43-5 




2.5x 10" 


2 


4.5x 10" 


3 


1.3x 10" 


6 


M8-1 




1.9x10" 


2 


2.2x10" 


3 


5.1x10" 


7 


M8-2 




1.3x10" 


2 


1.6x10" 


3 


5.1x10" 


7 


M8-3 




2. Ox 10" 


2 


2.6x10" 


3 


6.1 X 10" 


7 


M8-4 




2.0x10" 


2 


2.2x10" 


3 


5.5x10" 


7 


M8-5 




2.9x10" 


2 


2.6x10" 


3 


6.7x10" 


7 


M8-6 




1.6x10" 


2 


2.5x10" 


3 


5.5x10" 


7 


M16-1 




3.7x10" 


3 


1.0x10" 


3 


2.4x10" 


7 


M16-2 




2.0x10" 


2 


3.1x10" 


3 


4.6x10" 


7 


M17-1 




2.8x 10" 


-2 


2.6x 10" 


3 


7.7x 10" 


7 


Ml 7-2 




3.4x10" 


2 


3.5x10" 


3 


1.2x10" 


6 


Ml 7-3 




4.0x10" 


2 


3.1x10" 


3 


1.0x10" 


6 


M20-1 




6.3x10" 


2 


1.3x10" 


2 


2.2x10" 


6 


M20-2 




1.8x10" 


2 


3.0x10" 


3 


6.3x10" 


7 


M20-3 




7.6x10" 


3 


2.3x10" 


3 


3.6x10" 


7 


NGC 7635- 


-1 


6.5x10" 


3 


1.0x10" 


3 


2.7x10" 


7 


NGC 7635- 


-2 


1.1x10" 


2 


1.5x10" 


3 


4.1x10" 


7 


NGC 7635- 


-3 


9.1x10" 


3 


1.5x10" 


3 


3.4x10" 


7 



Derived from Fe/0 x 0/H 



Fe (Rodriguez 1999b ). In fact, the only features show- 
ing distinct behaviour in the spectrum of M20-1 are the 
[Feiii], [Fell] and [Niii] hues, whose intensities relative 
to H/3 are significantly higher than those measured for 
M20-2 and M20-3 (see Tables 1 and 2). Even the relative 
intensities of the [Fell] lines, whose values reflect the ef- 
fects of fluorescence in most [Fell] lines (Rodriguez 1996, 



1999a), depart in M20-1 from the line ratios measured 
for the other objects. For example, the intensity ratio be- 
tween [Fell] A4287, expected to be very sensitive to fluo- 
rescence, and [Fell] A8617, almost insensitive to pumping 
effects, is found to be lower in M20-1 than in any other 
position, suggesting that fluorescence effects are minimal 
in this position. Therefore, although [Fell] A8617 could in 
principle be affected by fluorescence in some objects, the 
high values derived for Fe+/H+ and Fe"^/Fe+"^ in M20-1 



(the value of this latter ratio being above those found for 
lower-excitation objects like M43 and NGC 7635), cannot 
be explained by resonance effects. The [Fell] A8617 line 
could only be measured in another position in this nebula 
(M20-2), where Fe+/Fe++ is found to be equally high, 
but since this line intensity is highly uncertain in M20- 
2, it is not clear whether these values of Fe+/Fe++ are 
characteristic of M20 or are peculiar to the jet sampled in 
M20-1 and imply a higher dust destruction efficiency in 
the lower ionization layers of the jet. In any case, a sim- 
ilar, although less extreme, Fe overabundance in M20-1 
can be found when the Fe abundances are calculated from 
just [Fc III] emiss ion an d model-predicted ionization frac- 
tions (Rodriguez 119961) . This overabundance of Fe atoms 
in the gas, which has also been found for several ty pical, 
not fully ionized HH objects (Beck-Winchatz et al. 1996 



and references therein), presumably reflects dust destruc- 
tion due to shock waves arising from the interaction with 
the ionized gas of the material in the jet, which is esti- 
mated to move at a velocity of ~ 370 km s~^ (Hester et 
al. |1999D F|. 



7.2. [Feiv] A2837 and Fe3+ abundance 

The only available measurement of an [Feiv] line in an 
Hii region is that performed by Rubin et al. ( 1997 ) for 
[Feiv] A2837 in M42. Their measured reddening-corrected 
intensity ratio, J(2837)//(H/3) ~ 10"^, along with two 
previous photoionization models for two regions in M42 
and the collision strengths calculated by Berrington & 
Pelan (|l99|, |l996| ) for Fe^+ , imply Fe/H = 4.6 x 10"'^ 
and 1.6 x 10^^ (with the difference between the two mod- 
els' results being due mostly to their different predicted 
temperatures). These abundance ratios are significantly 
below Fe/H = 3.0 x 10"^ the value Rubin et al. (|1997D 



consider consistent with [Fe ii] and [Fe ill] emission in M42 
and their models for this nebula. 

The discrepancy is somewhat reduced if the abun- 
dances derived here from [Fell] and [Felli] lines (and 
more recent atomic data for [Feiii]) in M42 are consid- 
ered: Fe/H = (0.8-1.8) x lO"*^. Although the weakness of 
the only measured [Fe iv] line, its sensitivity to tempera- 
ture and the extinction correction may account for part 
of the discrepancy, it could be entirely due to errors in 
the collision strengths for Fe^+. However, the more ex- 
tensive calculations of the collision strengths for this ion 



performed by Zhang & Pradhan (1997) lead to very simi- 
lar abundances. A 33-level model atom of Fe^+ based on 



the collision strengths of Zhang & Pradhan (1997) and the 
transition probabilities from Garstang (195S), for temper- 
atures in the range 8000-8500 K and densities of several 
thousands (these physical conditions being in agreement 
with those usually measured in different zones in M42 from 
[O III] and [CI ill] lines that should originate in the [Fe iv] 
emitting region: e.g. McCall 1979; Torres-Peimbert et al. 
1980| Esteban et al. |199S| ; Rodriguez |l9991j ) and the in- 



See also jittp:/ /eagle. Ia.asu.edu/hester/m20. html 
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tensity ratio measured for [Feiv] A2837 by Rubin et al. 
( |1997D , gives Fe^+ ~ 2 x 10"'^. This value is a fac- 
tor of 2-8 lower than the Fe'^+ abundances implied by 
the results presented here for M42 in Fig. ^ and Table 3: 
Fe'^^/H"'" — (4-15) X 10^^. The discrepancy is merely of a 
factor of - 2 if only M42 A-1 and M42 A-2, located only 



8'.'5 south of the position observed by Rubin et al. (1997), 
are considered. Given all the uncertainties related to the 
determination of Fe'^+/H+ from just one UV, weak [Fe iv] 
line measured in only one position whose optical spectra 
and [Fell] and [Feiii] line intensities are not available at 
the moment, and given the wide range covered by the val- 
ues derived here for Fe/H in M42 and their uncertainties, 
a discrepancy of a factor > 2 might not be significant. 

Taken at face value, the discrepancy could indicate a 
failure in the model predictions for the Fe ionization frac- 
tions and a corresponding uncertainty in the Fe abun- 
dances derived for high- ionization Hii regions. It has 
also been suggested that the lower Fe abundance im- 
plied by the [Fe iv] line could be real, reflecting the pres- 
ence of a gradient in the gaseo us abundance of this ele- 
ment (Bautista & Pradhan 1998 ), but the uncertainties re- 
lated to the determination of the Fe'^^ abundance are cur- 
rently too high to reach any definite conclusion. Further 
studies of [Feiii] and [Fell] lines at the position where 
[Feiv] A2837 has been observed and an assessment of the 
accuracy of the atomic data for Fe'^^ will be needed to 
reach a full understanding of [Fe iv] emission in H ii re- 
gions. 

8. Fe abundance, extinction and dust in H ii 
regions 

8.1. Fe abundance and dust destruction 

The high depletion factors measured for Fe in the inter- 
stellar medium and the high cosmic abundance of this ele- 
ment imply that Fe must contribute a substantial amount 
of mass to the dust grains (Sofia et al. 1994). Therefore, 
the degree of incorporation of Fe atoms to dust should 
give a good indication of the dust-to-gas ratio. Since the 
seven H II regions studied here have similar chemical abun- 
dances for those elements usually found in the gas phase 
(Rodriguez 1999b), it will be assumed that the Fe total 
abundance, i.e. in gas and dust, is the same for all these 
objects. The Fe abundance values derived here can then 
be compared to a reference Fe abundance to determine the 
fraction of Fe atoms present in the gas and depleted into 
dust for each area considered. Taking the solar Fe abun- 
dance as reference, the values derived for Fe/H or Fe/N 
imply that about 0.8-6% or 1-19% of the Fe atoms have 
been released to the gas phase. The corresponding log- 
arithmic depletion factors, shown in Fig. ^, fall between 
the values typically found for cool clouds in the Galactic 
disc, where log(Fe/H) - log(Fe/H)0 ~ -2.09 to -2.27, 
and those measured in the Galactic halo, whic h are in the 
range ~ -0.58 to -0.69 (Savage & Sembach |l996|) . The 
relatively high Fe abundances in the halo are similar to 



the values found for high- velocity clouds, and this differ- 
ence in the Fe depletion factors for cool disc and halo 
clouds (and the intermediate values found for warm disc 
clouds and clouds sampling both the disc and the halo) 
is interpreted as reflecting the action of the primary dust 
destruction agent: shock waves due to supernovae. 



The Fe abundances derived here imply that most of 
Fe remains in the dust grains inside Hii regions, but also 
imply the release to the gas of a significant amount of 
these atoms from the dust grains originally located in- 
side molecular clouds, where up to 100% of Fe is expected 
to be in the solid phase. Energetic photons, instead of 
shock waves, (excluding M20-1, see Sect. 7.1) seem to be 
responsible for this dust destruction effect, as suggested 
by the correlation of the Fe abundances with the degree 
of ionization in Fig. |3[ Two processes might be respon- 
sible for this effect: vaporization and optical erosion or 
photodesorption, both leading to the release of atoms or 
molecules from the grains after the absorption of photons. 
Optical erosion results from the ionization or excitation 
of a bonding electron of a molecule after the absorption 
of a photon; this mechanism is expected to destroy only 
the ice mantles on the grains, and the refractory grains 
should not be affected (Ostcrbrock 1989| ). [However, note 
that photodesorption cross-sections for metal atoms on 
dust grains are not known, and this process could be of 
importance.] Assuming then that sublimation is the pro- 
cess that releases Fe atoms to the gas, the grains contain- 
ing these atoms cannot be classical, i.e. characterized by 
a homogeneous composition and structure and with radii 
a ~ 0.1 fira, since these classical particles are in thermal 
equilibrium with the radiation field and have tempera- 
tures T ^ 100 K according to the calculations and the 
estimates obtained from the observed infrared emission, 
whereas refractory particles have much higher evapora- 
tion temperatures: T > 1000 K. However, if the grains 
are small enough, a < 0.03 fim, they can undergo sig- 
nificant temperature fluctuations after the absorption of 
one or a few energetic photons. Several studies on the 
temperature fluctuations undergone by dust grains in dif- 
ferent environments conclude that only very small grains, 
with a ~ 10 A, reach temperatures high enough to be de- 
stroyed. As an example, Guhathakurta & Draine (1989) 
estimate that graphite or silicate grains with a ~ 5 A lo- 
cated at 0.3 pc from a B3V star could reach temperatures 
T ^ 2000 K, would lose atoms through sublimation and 
have lifetimes of about 3 x 10^ years (whereas the ages of 
Hii regions are of the order of one million years). These 
small dust grains contain about one thousandth of the to- 
tal dust mass (Sellgren 1984) and their destruction cannot 
explain the relatively high Fe abundances measured for 
M42. Nevertheless, a further population of small grains 
could be located in larger composite grains. Interstellar 
grains, and in particular those from molecular clouds, are 
expected to consist of loose aggregates of smaller subunits 
formed by an accretion process that might fragment once 
inside the ionized gas or at the ionization front. 
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As an example, PAH-like clusters within larger amor- 
phous carbon grains could be easily released. Note that 
PAHs - polycyclic aromatic hydrocarbons, molecules of 
size ^ 10 A proposed to explain several emission bands 
observed in different environments such as the ionization 
fronts in Hil regions - are destroyed inside Hii regions 



(Bregman et al. |1989[ Giard et al. |l99j). PAHs or other 
kinds of aromatic molecules could aggregate and form 
conglomerates or amorphous carbon grains inside dense 
molecular clouds, and the abundant and chemically active 
Fe atoms might react with PAHs and accelerate the pro- 
cess. In fact, organometallic complexes could contain 5 to 
10% of metal atoms in the interstellar medium and, thus, 
these complexes would play an important role in the de- 
pletion of heavy elements and in dust fo rmati on processes 
in th e interstellar medium (Serra et al. 1992 ; Klotz et al. 



1995) 



8.2. Fe abundance and colour excess 

The Fe/0 abundance ratios are shown as a function of 
the colour excess implied by the relative intensities of the 
Hi lines in Fig. ^ where it can be seen that there is a 
loose correlation between the two quantities. The corre- 
lation of Fe/H and Fe/N with E{f3 — a) is also appar- 
ent, although somewha t less markedly than that of Fe/0. 
Previously, Rodriguez ( 1996| ) presented the correlation of 
Fe/0 with Tp (the effective optical depth for extinction at 
H/3), but the colour excess, which is an observed quantity, 
independent of the extinction law, is a more reliable and 
straightforward variable. 

The positions in M17 do not follow the correlation 
in Fig. ^, but this region is known to be heavily red- 
dened by foreground dust, both inters tellar and associated 
with the complex: Chini & Wargau ( 1998 ) estimate that 
E{B—V) ^ 1 mag is due to the diffuse interstellar medium 
in the line of sight to M17 (note that £'(/?-«) ~ E{B-V) 
when Rv — Av/E{B — F) ~ 3, the typical value for the 
diffuse interstellar medium) whereas the internal stars are 



further affected by Ay > 3 mag of local extinction. The 
positions studied in the other Hii regions might be only 
slightly affected by external extinction due to a selection 
effect, since the observations were centred on the bright- 
est positions of objects with high superficial brightness, 
but the difficulties inherent in disentangling the extinc- 
tion components due to external and internal dust, pre- 
clude an estimate of the effects of interstellar extinction 
(i.e. unrelated to the complex) in the correlation presented 
in Fig. ^. 

The extinction due to the interstellar medium along 
the line of sight to an Hii region can be estimated by 
studying the stars projected on the nebula and assuming 
that those whose extinction law is similar to that usually 
found for the diffuse interstellar medium should be located 
in front of the stellar formation regio n. This is the method 
followed by Chini & Wargau ( [1998D for M17, and i t has 
also been applied to M42, where Breger et al. ( 1981 ) find 
that E{B — V)~ 0.05 mag is due to interstellar extinc- 
tion. The distinction between extinction due to internal 
dust and that due to external, associated dust is more dif- 
ficult to make. However, taking into account the effects 
of scattering by dust internal to H ii regions, this internal 
dust mixed with the ionized gas can produce a maximum 
colour excess E{(3 — a) ~ 0.19 mag (Gomez Garrido & 
Miinch 1984). Therefore, the correlation in Fig. |^ relates 
the Fe abundance in the ionized gas - arising from dust 
destruction taking place in the region considered - with a 
colour excess necessarily due to foreground dust. 

Given a volume of gas and dust, the spatial separa- 
tion between the two components could be brought about 
through the effects of radiation pressure or stellar winds 
on the dust grains. However, coUisional interactions and 
Coulomb forces between ions and charged particles tend 
to hold dust and gas together in Hii regions. The effi- 
ciency of these processes depends on the size, structure 
and composition of the grains; the available calculations 
show that classi cal pa rticles are strongly coupled to the 
gas (Osterbrock 1989 ). S mall grains could be an excep- 
tion: Cardelli & Clayton (1988) estimate that grains with 
sizes a — 0.01 /im initially located at a distance r from the 
central star can move under the effects of radiation pres- 
sure by H^r/r > 0.2 in 10^ years. Since many Hii regions 
are known to have complex structures, with dense clumps 
or molecular columns inside or in front of the ionized gas 
(e.g. MB and NGC 7635, as shown by their detailed HST 
images), the overall effect of radiation pressure or stellar 
winds could well be to push a significant amount of dust 
outside the nebula into our line of sight. 

In conclusion, both correlations shown by the gaseous 
Fe abundance in H ii regions seem to require the presence 
of dust grains differing from the classical grains due to 
their smaller sizes or loose structure. These non-classical 
grains could contain 10-20% of the total Fe abundance, 
their destruction would be due to photons with energies 
above 35 eV and could take place after the fragmentation 
of an original composite grain. On the other hand, on a 
longer timescale, the surviving non-classical grains would 
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be pushed outside the ionized region through the action 
of radiation pressure or steUar winds, where they would 
be responsible for most of the colour excess found in the 
objects studied here. 

9. Conclusions 

The spectra of several positions inside seven Galactic H ii 
regions have been analysed to derive their Fe abundance. 
The relative intensities measured for the [Fe ill] lines have 
been compared with the values predicted by three com- 
binations of the available atomic data, and it has been 
shown that the collision strengths of Zhang (|1996D and 
the transition probabilities of Quinet ( 1996 ) lead to the 
best agreement between the observed and predicted line 
ratios. The Fe"'""'"/H"'" abundance ratios derived for all the 
positions studied with the above mentioned atomic data 
are 10-50% lower than t he va lues implied by the calcu- 
lations of Keenan et al. ( 1992 ), which are based on older 
atomic data and have been used in several previous studies 
of the Fe abundance in nebulae. 

The Fe+ZH"*" abundance ratios have been calculated 
using the [Fell] A8617 line, which is almost insensitive to 
fluorescence, and the cmissivities of Bautista & Pradhan 
( |1996| ). The derived Fe+ abundances may be affected by 
uncertainties in the measured intensities and in the colli- 
sion strengths, but the effects of these uncertainties in the 
Fe abundances should be small, since the contribution of 
Fe"*" to the total abundance is very low for most posi- 
tions. The Fe abundances are consequently derived from 
the values of Fe"'""'"/H+ -I- Fe"''/H+ and ionization correc- 
tion factors for the contribution of Fe^"*" obtained from 
available grids of photoionization models. 

The derived Fe/0 abundance ratios are 3 to 40 times 
lower than the solar value, implying that most of the Fe 
atoms remain in dust grains inside the ionized volume. 
These underabundancc factors are, however, significantly 
lower than those found for the undisturbed interstellar 
medium, where more than 99% of Fe is depleted into dust 
(Savage & Sembach 1996 ). The correlation between the 
calculated Fe abundances and the degree of ionization 
explains this difference by relating the energetic photons 
in Hii regions with the release of Fe atoms from dust. 
Further agents can be responsible for dust destruction in 
H II regions, and a remarkable example is provided by po- 
sition M20-1, on the tip of the optical jet HH 399, where 
the Fe gaseous abundance is significantly higher than in 
the other areas, presumably indicating dust destruction 
due to the shock waves related to the jet. 

The inferred dust destruction by photons, together 
with the correlation of the Fe abundance with the colour 
excess, can be explained by invoking the presence of 
loosely-bound dust grains or mantles, containing 10 to 
20% of the Fe abundance, whose fragmentation would 
produce grains small enough to be destroyed by sublima- 
tion processes, whereas the remnants of these small grains 
would be swept out of the ionized region by radiation pres- 
sure or stellar winds. 



The present results rely on the estimated contribution 
of Fe'^^ to the total Fe abundance. The uncertainties re- 
lated to this estimate, strengthened by the current difficul- 
ties in explaining [Feiv] emission in M42, may cast some 
doubts on the validity of the results for the high ionization 
objects. Certainly, this issue should be resolved before the 
correlation between the Fe abundances and the degree of 
ionization can be considered as firmly established, thus 
favouring the proposed interpretation. If such were not 
the case, and the predicted contribution of Fe"^+ to the 
total abundance were found to be grossly in error, the 
new results would surely lead to alternative scenarios for 
the evolution of dust in H ii regions. 
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